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Abstract 

The BaBar detector, at the PEP-11 asymmetric B Factory at SLAC collected a sample of 32 
million BB pairs whilst operating at energies near the T[AS) resonance between October 1999 and 
May 2001. An study of time-dependent CP- violating asymmetries in events where one neutral B 
meson is fully reconstructed in a final state containing charmonium produced the measurement 
sin 2/3 = 0.59 it 0.14 (stat) it 0.05 (syst), which constitutes an observation of CP violation in the B^ 
meson system at the Aa level. Also presented are preliminary results from a study of CP violation 
in the decays B^ —>■ 7r"'"7r~ and B^ K~^Tr~. 
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1 Introduction 



CP violation was first observed in the decays of Kj^ mesons in 1964 It was 37 years before this 
phenomenon was observed in another system, that of the neutral B mesons. The effect arises from 
the existence of an irremovable, CP-violating phase in the three-generation CKM quark-mixing 
matrix [Q]. The system has an important advantage over the system in that the measure- 
ments of CP-violating asymmetries provide a direct test of the Standard Model of electroweak 
interactions, free of corrections from strong interactions, which arise when trying to interpret the 
results from decays. The primary goal of BaBar is to over-constrain the Unitarity Triangle 
through multiple, independent measurements of its sides and angles. In this analysis, the angle (5 
is probed through a measurement of sin 2/? using b — > ccs decays. 

2 PEP-II and BaBar 

The PEP-II B factory at SLAC comprises a pair of storage rings producing asymmetric e~^e~ 
collisions (9 GeV e~, 3.1 GeV e"^) at a center-of-mass energy corresponding to the mass of the 
T(4S) resonance (10.58 GeV). The T(4S) decays almost exclusively to B~^B~ or coherent B^B^ 
pairs. The primary goal of the BaBar detector for this analysis is to measure the time difference 
between the two meson decays, At. The asymmetric nature of the collisions gives a boost of 
/37 = 0.55 which yields an average spatial vertex separation of ~ 250 ^m. 

A detailed description of the BaBar detector can be found in Q. The volume inside the 1.5T 
superconducting solenoid consists of a five layer silicon vertex detector (SVT), a drift chamber 
(DCH), a quartz Cerenkov detector (DIRC) and a CsI(Tl) crystal electromagnetic calorimeter 
(EMC). The instrumented flux return (IFR) outside the magnet comprises alternate layers of iron 
and resistive plate counters (RPCs). 

3 Exclusive B reconstruction 

A sample of neutral B mesons, Bcp have been fully reconstructed in their decays to final states 
of known CP content: J/i^K'^, ij{2S)K'^, J/ipKl, XciK^ and J/i}K*^{K*^ K^t^^). In addition 
there are also samples of B decays to final states of definite flavour (iJ^av): B^ D^*^~tt~^, 
D^*^-p+, D^*^'at and J /il}K*^{K*^ K+-K-) as weh as charged B decays: B~ L»(*)°7r-, 
J/^K~ and ip{2S)K^ (charge conjugate modes are implied throughout this paper). All selections 
have been optimized to give maximum sensitivity to the final measurement. Particle identification, 
mass (or mass difference) and vertex constraints are used wherever applicable. The resulting 
signal yield for each mode is identified by using the kinematical variables /S.E = E*^ — E^^^^ 

and ruES = \J ^leam ~ P*B where E'^, p*^ are the center-of-mass energy and momentum of the 
reconstructed B and E^^^^^ is the beam energy in the center-of-mass. In the case that an event has 
more than one B candidate, only the one with the smallest | A£^| is retained. For each mode a signal 
region is defined as ±3(7 about (5.279,0) in the tues, plane. The niEs resolution is f« 3MeV/c^, 
dominated by the spread of the beam energy. The A.E resolution is mode dependent and varies 
from about 10-33 MeV. Figure |^ shows the mEs distributions for Bcp candidates containing a 
and the A£^ distribution for candidates containing a K^. The number of tagged events and the 
signal purities, determined from fits to the mES (all modes except K*^) or /S.E {K^ mode) 
distributions in data or from Monte Carlo simulation {K*^ mode) are shown in table |l[ 
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Figure 1: a) Distribution of ruEs for Bcp candidates having a Kg in the final state; b) distribution 
of AE for J/ipKf^ candidates. 



4 Flavour Tagging 

Flavour tagging information is extracted from the other (partially reconstructed) B in the event, 
Stag- The coherent production of the B^B^ pair ensures that the flavour of Bcp is exactly opposite 
to that of Stag at the time when Stag decays, ttag- Each event is assigned to one of four hierarchical, 
mutually exclusive tagging categories or excluded from further analysis. The Lepton and Kaon 
categories contain events with high momentum leptons from semileptonic B decays or with kaons 
whose charge is correlated with the flavour of the decaying h quark (e.g. a positive lepton or kaon 
yields a B^ tag). The NTl and NT2 cate gories are based on a neural network algorithm whose 
tagging power arises primarily from soft pious from D*^ decays and from recovering unidentified 
isolated primary leptons. The Sflav sample is used to measure the tagging performance along with 
the Bcp events. The figure of merit used \s Qi = €i{l — 2wi)'^ where ei and wi are the efficiency 
and mistag fraction for category i. The statistical error on sin 2/? is proportional to where 
Q = Y1, Qi- The efficiencies and mistag fractions for the four tagging categories are shown in table |2|. 



Table 1: Number of tagged events, signal purity and result of fitting for CP asymmetries in the 
full CP sample and in various subsamples, as well as in the -Bflav and charged B control samples. 
Errors are statistical only. 



Sample 


^tag 


Purity (%) 


sin 2/3 


J/^ljK^„ij{2S)KlxciK'', 
J/i)K*^,K*^ i^V 


480 
273 
50 


96 
51 
74 


0.56 ± 0.15 
0.70 ± 0.34 
0.82 ± 1.00 


Full CP sample 


803 


80 


0.59 ± 0.14 


J/vK'': ti2S)K'^, AciA"!* oiib 


' ('//■ 


= -1) 




J/V^ir^ (ir^ ^ TT+TT-) 
J/T/^fSfO ^ 7tOtT°) 
?/;(25^)i^:° {K^g TT+TT") 
Xci-f^s [J^s ~^ ) 


316 
64 
67 

OO 


98 
94 
98 
y / 


0.45 ± 0.18 
0.70 ± 0.50 
0.47 ± 0.42 

4- 0-55 


Lepton tags 


74 


100 


0.54 ± 0.29 


Kaon tags 


271 


98 


0.59 ± 0.20 


NTl tags 


46 


97 


0.67 ± 0.45 


NT2 tags 


89 


95 


0.10 ± 0.74 


tags 


234 


98 


0.50 ± 0.22 


tags 


246 


97 


0.61 ± 0.22 


Bflav non-CP sample 


7591 


86 


0.02 ± 0.04 


Charged B non-CP sample 


6814 


86 


0.03 ± 0.04 



5 measurement and resolution 

The time difference between the two B decays, At = tcp — Uag is determined by first measuring 
their spatial separation Az = zcp — ztag- This is corrected on an event-by-event basis for the 
direction of the B with respect to the z direction in the T{4S) frame, zcp is determined from 
the charged tracks that constitute the Bcp candidate. All other tracks in the event are fitted to a 
common vertex in order to calculate the i^tag decay position, Ztag- Tracks from photon conversions 
are removed. Pairs of tracks compatible with the decay of a long lived Kg or A are replaced by their 
parent neutral pseudotrack. The bias in the forward z direction due to charm decays is reduced 
by removing the track with the largest contribution to the vertex x^, if above 6 and iterating the 
fit until either no track fulfills this condition or fewer than two tracks remain. Knowledge of the 
beam spot location and beam direction is incorporated through the addition of a pseudotrack to 
the tagging vertex, computed from the (-Bflav) vertex and three-momentum, the beam spot 
(with a vertical size of 10 /xm) and the T{AS) momentum. The total Az reconstruction efficiency 
is 97%. For 99% of the reconstructed vertices the r.m.s. Az resolution is 180 /um, dominated by 
the Stag vertex. An accepted candidate must have a converged fit for the Bcp and i?tag vertices, 
an error of less than 400 /xm on Az and a measured |Ai| < 20 ps. 

The At resolution function for signal events is represented as a sum of three Gaussian dis- 
tributions. All offsets are modelled to be proportional to the event-by-event error, (JAt, which is 
correlated with the weight that the daughters of long-lived charm particles have in the tag vertex 



Table 2: Efficiencies and average mistag fractions Wi extracted for each tagging category i from a 
maximum likelihood fit to the time distribution for the fully reconstructed B sample (Sep + i?flav)- 
Uncertainties are statistical only. 



Category 


6(%) 


w{%) 


Q{%) 


Lepton 


10.9 ±0.3 


8.9 ± 1.3 


7.4 ±0.5 


Kaon 


35.8 ±0.5 


17.6 ±1.0 


15.0 ±0.9 


NTl 


7.8 ±0.3 


22.0 ±2.1 


2.5 ±0.4 


NT2 


13.8 ±0.3 


35.1 ±1.9 


1.2 ±0.3 


All 


68.4 ±0.7 




26.1 ± 1.2 



reconstruction. The 'core' and 'tail' Gaussians have widths scaled by the event-by-event measure- 
ment error derived from the vertex fits. A separate offset for the core distribution is allowed for each 
tagging category to account for small shifts caused by inclusion of residual charm decay products 
in the tag vertex. The third Gaussian has a fixed width of 8 ps and accounts for fewer than 1% of 
events with incorrectly reconstructed vertices. Identical resolution function parameters are used for 
all modes, since the At resolution is dominated by the -Btag vertex precision. Separate resolution 
function parameters have been used for data collected in 1999-2000 and 2001, due to the significant 
improvement in the SVT alignment. 

6 Measuring sin 2/3 

6.1 Improvements to the analysis 

There are several significant changes in this analysis |Q] relative to the first BaBar sin 2/3 publica- 
tion 1^]. The Bcp modes Xcil^s J/tpK*^{K*^ — > Kgir^) have been added. Improvements in 
track and Kg reconstruction efficiency in 2001 data produce a 30% increase in the yields per 
luminosity unit. Better alignment of the tracking systems in 2001 data and improvements in the tag 
vertex reconstruction algorithm have increased the sensitivity of the measurement by an additional 
10%. The purity of the sample has been increased by a reoptimization of the J /ipK^ selection. 
In total, the statistical power of the analysis is almost doubled with respect to that of Ref. ||5|. 
The final Bcp sample contains about 640 events, the -Bflav sample has 7591 fully reconstructed B^ 
events and the charged B sample has 6814 fully reconstructed B^ events. 

6.2 The sin 2/? fit 

The fit for sin 2/3 is based on the following framework. Each event in the Bcp sample is examined 
for evidence that the Stag decayed as a B^ or B^. The decay distributions as a function of time for 
events with either type of tag can be expressed in terms of a complex parameter A that depends 
both on B^B^ mixing and the amplitudes describing B^ and B^ decay to a common final state 
/ @. The distribution f+(f_) of the decay rate when the tag is a B^{B^) is 



^ 2rso(l + |A|2) [ 2 
1 - lAp 

lb TmAsin (Am^oAt) =p cos (Am^oAt) 

where t^o is the hfetime and AM^o is the oscillation frequency of B^B^ mixing. The first 
oscillatory term is due to interference between direct decay and decay after mixing. A difference 
between the B^ and B^ At distributions or a At asymmetry for either flavour tag is evidence for 
CP violation. 

In the Standard Model, A = r/je~^*^ for h — > ccs decays containing charmonium where rjf is the 
CP eigenvalue of the state / and (3 = aigl—VcdVcb/ ~ ^tdVti^ is an angle of the Unitarity Triangle. 
Thus, the time-dependent CP-violating asymmetry is 

_ f+(At)-f_(At) 
^^'^^^ = f+(At) + f_(At) 

= — r/j sin2/?sin (Am^o At), 

where rjf = —1 for J/il)K^, i/j{2S)Kg and XciKg and for J/tpK^. Due to the presence of even 
(L=0,2) and odd (L=l) orbital angular momenta in the J/i/jK*^{K*^ — >■ Kgn^) system there can be 
CP-even and CP-odd contributions to the decay rate. When the angular information in the decay 
is ignored, the measured CP asymmetry in J/t/jK*^ is reduced by a dilution factor D± = 1 — 2R±, 
where R± is the fraction of the L=l component. It has been measured as R± = (16 it 3.5)% which, 
after acceptance corrections leads to an effective T]f = 0.65 it 0.07 for this mode. 

The Bcp and -Bflav samples are used together in the unbinned maximum likelihood fit for the 
extraction of sin 2/3. A total of 45 parameters are varied in the fit, including sin 2/3 (1), the average 
mistag fraction w and the difference Aw between B^ and B^ mistags for each tagging category (8), 
parameters for the signal At resolution (16) and parameters for background time dependence (9), 
At resolution (3) and mistag fractions (8). The determination of the mistag fractions and signal At 
resolution function is dominated by the large -Bfiav sample. Background parameters are governed 
by events with tues < 5.27GeV/c^. As a result, the largest correlation between sin 2/3 and any 
linear combination of the other free parameters is only 0.13. The B lifetime and mixing frequency 
are fixed at r^o = 1.548 ps and AMgo = 0.472 ?ips~^ respectively ||7|. The value of sin 2/3 and the 
asymmetry in the At distribution were hidden using a blind analysis method (as described in ||^]) 
until the event selection was optimized and all other aspects of the present analysis were complete. 

Figure ^ shows the At distributions and Acp as a function of At overlaid with the likelihood 
fit results for the rjf = —1 candidates. The probability of obtaining a lower likelihood value than 
the one from the fit, evaluated using a Monte Carlo technique, is 27%. The simultaneous fit to the 
full sample yields 

sin 2/3 = 0.59 ± 0.14 (stat) ± 0.05(syst). 



6.3 Systematic errors and consistency checks 

The systematic error is dominated by the parametrization of the At resolution function (0.03) - 
due in part to residual uncertainties in SVT alignment, possible differences in the mistag fractions 
between the Bcp and -Bfjav samples (0.03) and uncertainties in the level, composition and CP 
asymmetry of the background in the selected CP events (0.02). 
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Figure 2: a) Number of ry^^ = — 1 candidates in the signal region a) with a tag, A^^o and b) with 
a tag, N^o and c) the asymmetry {N^o — N^o)/{Ngi) + N^o), as functions of At. The solid 
curves represent the result of the combined fit to all selected CP events and the shaded regions 
represent the background contributions. Figures d)-f) contain the corresponding information for 
the r]f = +1 mode. The likelihood is normalized to the total number of B^ and B^ tags. The value 
of sin 2/3 is independent of the individual normalizations and therefore of the difference between 
the number of B^ and B^ tags. 



The large sample of reconstructed events allows a number of consistency checks, including 
separation of the data by decay mode, tagging category and -Btag flavour. The results of fits to 
these subsamples and to the samples of non-CP decay modes (where no statistically significant 
asymmetry is found) are shown in table ||. 

7 Results from tt+tt , iC+Tr decays 

A search for CP-violating asymmetries in decays of P" to two light mesons has been carried out 
on a sample of 33 million BB pairs [^. Here one hopes to measure sin 2a although the analysis 
is complicated by the possibility of penguin pollution leading to diagrams with different weak and 
strong phases contributing to the same final state. In this case, one has to allow for direct CP 
violation and in general |A| 7^ 1. The fit is then complicated by the presence of a cosine term as 
well as the coefficient of the sine term containing sin2aeff) where a^s depends on the magnitudes 
and strong phases of the tree and penguin amplitudes. The selected signal B sample consists of 
65j^|^ vrvr, 217 it 18 Ktt and 4.3^*^43 events. The results for CP-violating asymmetries are 
summarized in table ^ Here S^tt and C^n are respectively the coefficients of the sine and cosine 
terms in the expression for f-i-(At) and 

Ak = ^^"^^ ~ 

Figure |3| shows the At distributions and the asymmetry Amri^t) = (N^o {At)—N^o i^t))/ {N^o (At)+ 
N^o{At)) for tagged events which are enhanced in signal tttt decays. 



Table 3: Central values and 90% confidence level intervals for S-^t^, C-j^n and Ak- 



Parameter 


Central Value 


90% C.L. Interval 


c 


o.o3l[J:i±o.ii 

-0.25lgf ± 0.14 
-0.07 ±0.08 ±0.02 


[-0.89,±0.85] 
[-1.0,±0.47] 
[-0.21,±0.07] 



8 Conclusions and outlook 

The measurement of sin 2/? presented here establishes CP violation in the B^ meson system at 
the 4.1(7 level, 37 years after its discovery in the Kaon system. The probability of obtaining this 
value of sin 2/3 or higher in the absence of CP violation is 3 x 10~^. This direct measurement is 
consistent with the range implied by measurements and theoretical estimates of the magnitudes of 
CKM matrix elements By the summer of 2002, with a data sample of more than 100 million 
BB pairs a measurement of sin 2/3 with a precision of less than 0.1 will be possible. 

In addition the search for CP-violating asymmetries in the decays B^ — > vr+vr" and B^ — > K~^tt~ 
looks promising and with a similar amount of data should yield errors of ~ 0.3 on the B^ tt^tt" 
asymmetries. 




Ar (ps) 

Figure 3: Distributions of At for events enhanced in signal vrvr decays. Figures (a) and (b) show 
events (points with errors) with -Btag = or B^. Solid curves represent projections of the max- 
imum likelihood fit, dashed curves represent the sum of qq and Kir background events and the 
shaded region represents the contribution from signal vrvr events. Figure (c) shows ATmiAt) for 
data (points with errors) as well as fit projections for signal and background events (solid curve) 
and signal events only (dashed curve). 
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